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kThis  report  describes  several  investigations  performed  towards  improving 
our  ability  to  communicate  through  the  ionospheric  medium.  These  studies 
covered  two  distinct  and  separate  areas.  The  first  was  concerned  with 
improvements  in  ionospheric  forecasting  in  the  polar  region  through  the 
development  of  empirical  models  for  electron  density  profiles  and  iono- 
spheric parameters.  The  second  consisted  of  the  study  of  various  tech- 
niques for  the  establishment  of  high  frequency  communication  circuit  a 
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1.  The  final  report  on  "Ionospheric  Analysis  and  Ionospheric  Modelling", 
Contract  No.  F19628 -74-C-0070  consists  of  four  sections.  In  the 
modelling  effort  of  foF2,  daily  hourly  observations  of  six  years  of  the 
1958/68  time  span  have  been  incorporated  for  adjusting  the  ITS  model  to 
the  polar  region.  No  significant  improvement  over  the  ITS  model  is 
claimed  to  have  been  noted  after  testing.  In  ray-tracing  studies  for 
parabolic  and  quasi -parabolic  electron-density  height  variations  and 
horizontal  gradients  a transmitter -coordinate-system  was  introduced  to 
facilitate  the  assessment  of  the  effects  of  longitudinal  and  transverse 
gradients  on  propagation  characteristics.  An  iterative  scheme  was 
developed  to  place  a nominal  electron  density  profile  near  ray  reflec- 
tion. A functional  representation  of  effective  tilt  in  a parabolic  iono- 
sphere was  deduced.  The  three -dimens ional -model  effort  dealt  with  the 
synthesis  of  electron  density  height  profiles  by  means  of  a three-layer 
ionosphere  considering  monotonically  increasing  electron  density  profiles 
as  well  as  those  possessing  an  ionization  valley.  Procedures  for  proper 
matching  of  profile  segments  of  contiguous  layers  were  developed.  The 
coordinate  conversion  technique  for  OTH  backscatter  radar  applications, 
which  is  to  appear  in  a separate  report,  is  summarized.  A given  iono- 
sphere is  used  to  generate  the  leading  edge  of  a backscatter  ionogram. 

The  latter  is  compared  with  a given  backscatter  ionogram.  Range- 
gradient  adjustment  factors  are  applied  to  achieve  mutual  agreement. 

Once  agreement  is  reached,  predicted  radio  frequency  paths  are  simulated 
and  used  to  establish  the  relationship  between  ground  range  and  group 
path.  Eventually,  this  relationship  is  expected  to  yield  estimates  of 
ground  range  for  an  OTH-radar  target. 


2.  The  above  work  is  of  value  since  it  provides  computational  proce- 
dures for  utilizing  and  correcting  for  ionospheric  radio  propagation 
effects.  These  procedures  will  be  used  by  USAF  for  further  refinement 
of  the  ionosphere  model  to  more  realistically  predict  long-range  duct- 
ing in  ionospheric  channels  and  to  better  estimate  ground  range  to 
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I.  MODELING  foF2  IN  THE  POLAR  REGION 


In  this  section,  we  will  present  some  of  our  results  from  a study  to 
improve  our  knowledge  of  the  F2  region  of  the  ionosphere  in  the  polar 
region  using  an  empirical  model.  In  a previous  report^  \ we  described 
some  attempts  at  finding  more  accurate  models  of  foF2.  The  motivation 
and  method  of  approach  to  modeling  foF2  remain  the  same  as  in  our  earlier 
work. 

(2) 

The  ITS  monthly  median  model  for  foF2  has  been  shown  to  repre- 
sent the  physical  behavior  of  the  ionosphere  quite  well  in  the  midlatitude 
regions  of  the  globe.  We  have,  therefore,  chosen  to  use  this  model  as  a 
first  approximation  for  the  determination  of  foF2  at  high  latitudes.  Correc 
tions  to  the  model  are  then  applied  as  functions  of  the  relevant  spatial  and 
geophysical  variables.  While  the  ITS  model  is  only  based  on  median 
vertical  incidence  ionosonde  data  for  IGY195  8,  we  have  used  daily-hourly 
observations  for  1958,  I960,  1962,  1964,  1966  and  1968  in  the  determi- 
nation of  correction  terms  for  the  polar  region.  The  use  of  the  daily  data 
allowed  for  the  inclusion  of  the  planetary  magnetic  activity  index,  K 

P 

dependence  of  foF2  in  the  model.  Our  earlier  modeling  of  foF2  was  based 
primarily  on  the  1958  and  1964  data  and  the  resulting  changes  in  the  model 
reflect  in  large  part  the  use  of  available  data  from  other  years. 

The  ITS  monthly  median  model  represents  the  median  behavior  of 
foF2  in  terms  of  a spherical  harmonic  expansion  of  the  space  and  time 
variables.  The  representation  is  of  the  form 


N 

= ^ (T)  G^  (0,  X)  where  0 and 

k=0 


foF2 


1 


A are  the  geographic  longitude  and  latitude  and  T is  the  universal  time. 
Here  the  diurnal  variation  has  been  separated  from  the  spatial  variation 
and  the  cutoff  parameter  of  the  series,  N is  equal  to  75.  The  time  de- 
pendent functions  are  expanded  in  a Fourier  series  containing  six 
harmonics: 


Dk(T) 


+ Y (A.(k)  cos  jT  + 
J 

j = l 


B.(k)  sin  jT). 


This  results  in  a function 

N (k) 

foF2  = ^ Aq'  Gk(0,A)  + 
k=0 

Y T [A.(k)  cos  (jT)  Gk(0,\)  + B.(k)sin(jT)Gk(0,\)] 
k=0  j=l 

where  tie  G(0,  A)  are  linear  combinations  of  surface  spherical  harmonics. 
The  explicit  arguments  of  the  G functions  were  chosen  to  be  the  modified 
magnetic  dip  angle  x,  A and  0.  The  modified  magnetic  dip  angle  is  a 
function  of  the  latitude  and  magnetic  dip  angle  I,  given  by 


r 


B , B , B are  the  north,  east  and  vertical  components  of  the  earth 
x y z 

magnetic  field  vector.  The  modified  magnetic  dip  angle,  x,  is  then 


The  988  coefficients  were  determined  using  the  hourly  observations  of 
foF2  that  were  available  from  a network  of  vertical-incidence  ionosondes 
for  the  IGY  1958.  The  use  of  additional  data  from  1964  allowed  for  the 
determination  of  the  mean  sunspot  number  dependence  of  the  coefficients. 

The  stations  in  the  polar  region  for  which  foF2  data  was  available  for 
the  even  years  between  I960  and  1968  are  listed  in  Table  1-1.  Approx- 
imately 100,000  hourly  observations  of  foF2  were  available  for  each  of 
the  years.  While  the  stations  do  not  give  uniform  coverage  over  the 
polar  region,  they  can  be  used  to  explore  some  of  the  spatial  and  temporal 
features  of  the  F2  region  if  one  does  not  use  them  directly  with  simple 
regression  techniques. 

We  have  assumed  that  some  of  the  corrections  which  are  necessary 
for  the  polar  region  are  strongly  coupled  to  solar  phenomena  and,  there- 
fore, separate  all  observations  into  the  two  categories  of  solar  day  where 
\ (solar  zenith  angle)  £ 94.6°  and  solar  night,  \ > 94.6°.  In  addition  to  the 
foF2,  we  have  made  use  of  a data  base  of  values  of  the  three-hour  planetary- 

' 

magnetic-activity  index,  for  all  of  the  years  under  study. 

The  polar  ionosphere  can  be  described  in  terms  of  dynamic  effect 
which  occur  relative  to  several  spatial  boundaries.  The  most  notable  ones 
being  the  trough  and  auroral  oval.  In  our  previous  work,  we  assumed  that 


s 
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TABLE  1-1 


Stations  for 

Which  foF2 

Data  is  Available 

for  I960,  1962,  1964, 

1966,  1968 

Stations 

X (Degrees) 

X 

m (Degrees) 

0 (Degrees) 

Gorky 

156 

56.  2 

51.6 

44.  3 

Julensruk 

055 

54.6 

62.  5 

13.4 

Kiruna 

167 

67.  8 

64.  3 

20.  5 

Leningrad 

160 

60.0 

55.  8 

30.  3 

Lindau 

050 

51.4 

48.  3 

10.  1 

Lycksele 

164 

64.  6 

61.2 

18.  8 

Miedzeszyn 

152 

52.  2 

48.  2 

21.  2 

Moscow 

155 

55.5 

51.0 

37.  3 

Murmansk 

168 

68.  9 

64.  6 

33.  0 

Nurmijarvi 

159 

60.5 

56.6 

24.  6 

Providenya 

664 

64.4 

60.  3 

18.  5 

Sakekhard 

266 

66.  6 

61.  8 

66 . 6 

Slough 

051 

51.5 

49.  9 

359.4 

oodankyla 

1 66 

67.4 

63.4 

26 . 6 

Uppsala 

158 

59.  8 

56.  4 

17.  6 

Yakutsk 

462 

62.  0 

56.  4 

129.  8 

7 


the  boundaries  were  well  established  from  earlier  studies  and  examined 
the  observations  relative  to  these  specified  boundaries.  It  was  apparent 
from  our  data  that  the  spatial  position  and  magnitude  of  any  polar  region 
corrections  to  foF2  were  strongly  correlated  and  varied  widely  from  one 
year  to  another.  For  this  reason,  we  have  used  each  set  of  station-month 
data  to  determine  all  the  properties  of  any  feature  of  the  ionosphere  if  it 
exists  in  the  data.  In  addition  to  die  geographic  coordinates,  universal 
time,  and  sunspot  number  used  by  the  ITS  model  we  have  searched  for 
corrections  which  may  also  be  functions  of  the  corrected  geomagnetic 
space  and  time  coordinates  and  the  index.  As  a first  attempt  to  assess 
the  accuracy  of  the  median  model,  we  have  used  each  of  the  six  years 
foF2  data  to  calculate  the  RMS  residual^RMS  given  by: 


RMS 


, N 

N l 


, ITS 
foF2. 

1 


foF2. 


obs 


i=l 


. .-.ITS 
foF2. 
l 


where  the  index  i denotes  an  individual  measurement  for  a specific  universal 
time,  9,X,  sunspot  number  and  K^.  We  first  calculate  ^ IvIS  eac^  rnont^ 
and  for  solar  day  and  night  separately  averaging  over  all  station  positions, 
time  and  daily  geophysical  parameters.  The  results  are  shown  in  the  twelve 
curves  or  Figures  1-1  and  1-2. 

Although  the  median  model  does  contain  a dependence  on  solar  activity, 
the  data  presented  in  this  form  does  indicate  both  random  variations  and 
systematic  deficiencies  of  the  model.  The  solar  day  data  1966  appears  to 
give  a good  indication  of  the  minimum  random  RMS  error  of  the  model  of 


8 


k 


approximately  15%.  This  also  appears  to  be  of  the  same  order  in  1962 
and  19f>8,  but  there  is  some  evidence  of  systematic  effects  occurring 
during  the  spring  and  fall  equinox.  This  effect  is  much  more  pronounced 
during  the  other  three  years.  The  difference  between  the  solar  night  (SN) 
and  solar  day  (SD)  data  varies  greatly  over  the  years.  In  1964,  the  SN 
residual  is  almost  a constant  scale  factor  times  the  SD  residual  for  each 
month  of  the  year  while  the  SN  residuals  of  1968  exhibit  a strong  seasonal 
dependence  when  none  was  present  in  the  SD  residuals.  This  strong 
seasonal  nighttime  effect  is  present  in  all  the  years  with  the  ratio 

A™ » (winter)  / inw„(summer)  being  a factor  of  2 for  five  of  the  six  years. 

KMS  RMS 

It  should  be  emphasized  that  although  these  residuals  are  only  determined 
for  the  set  of  polar  region  stations,  the  coverage  of  the  region  may  still  be 
biased  with  respect  to  individual  dynamic  features  of  the  F2  region  of  the 
ionosphere. 

In  order  to  determine  the  corrections  to  the  ITS  model  for  foF2,  we 
have  divided  the  data  into  intervals  of  values  of  the  magnetic  activity  index 
and  the  corrected  geomagnetic  latitude.  Reasonable  statistics  were  obtained 
by  using  4°  intervals  of  geomagnetic  latitude  beginning  at  46°  North.  The 
1958  data  extended  to  74°  N while  all  the  other  years  only  gave  adequate 
coverage  to  62°  N.  The  data  was  divided  into  seven  intervals  of  magnetic 
activity  according  to  the  following  graph. 


* 


. ii 


For  convenience,  we  present  the  data  in  only  three  categories  of 
magnetic  activity: 

Low  K S 1.3 

P 

Medium  1 . 3 < K S 3,3 

P 

High  3.  3 < K . 

* P 

For  each  grouping  of  geomagnetic  latitude  and  magnetic  activity,  we 
calculate  the  average  percent  residual: 


A(X  ,K  ) 
m p 


-ill 


foF2. 


ITS 


foF2. 


obs 


i=l 


foF2J 


ITS 


where  the  summation  is  performed  for  all  stations  and  time.  The  results 

are  presented  as  a function  of  the  month  number  in  Figures  1-3  through  1-14 

for  each  of  the  years  in  the  study.  A strong  feature  of  almost  all  the  data 

is  the  similar  K dependence  of  the  average  residual  for  each  separate 
P 

latitudinal  interval.  The  correction  is  negative  for  low  magnetic  activity 
and  positive  for  high  magnetic  activity.  If  the  three  groups  are  averaged, 
the  results  are  consistent  with  the  median  model  which  was  formulated 
without  any  dependence  on  magnetic  activity. 

The  solar  day  and  night  residuals  for  1958  do  not  exhibit  any  noticeable 
dependence  upon  geomagnetic  latitude.  There  is  also  very  little  difference 
between  the  day  and  nighttime  curves.  This  is  consistent  with  earlier  work 
in  which  very  little  evidence  for  the  trough  and  auroral  oval  was  observed 
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AVERAGE  PERCENT  RESIDUAL:  SOLAR  DAY  1958 
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for  1958.  The  daytime  curves  for  the  other  years  also  show  very  little 
variation  as  a function  of  geomagnetic  latitude.  Since  this  data  is  confined 


to  thi-  ri-uion  of  X i (>L  , there  is  no  evidence  in  A duo  to  the  effect  of  the 
" m 


auroral  oval. 


In  the  curves  for  the  solar  night  data,  there  is  evidence  of  the  existence 
of  the  high-latitude  trough.  Although  it  is  not  present  in  all  years,  it  does 
show  a strong  seasonal  dependence  when  it  exists.  This  is  most  clearly 
seen  in  the  data  for  1968  where  the  ITS  model  predicts  a value  of  foF2 
which  is  approximately  25%  too  large  in  the  winter  months.  The  latitudinal 
dependence  of  the  trough  correction  is  also  visible  in  the  data,  but  appears 
to  be  only  slightly  correlated  with  magnetic  activity.  We  have  used  this 
data  to  estimate  the  position,  width  and  amplitude  of  the  trough  when  it  is 


observed.  For  a given  group,  year  and  month  the  average  residual  was 
examined  as  a function  of  X . For  each  month  that  shows  a peak  in  A(X  ), 


m 


m 


the  center  position,  the  full  width  at  one  quarter  maximum,  T and  the 


maximum  amplitude,  A was  estimated.  The  results  are  presented  in 

m 

Table  I-  2 for  the  four  years  I960,  1962,  1966  and  1968.  The  entries  in 
each  month  for  each  of  the  years  are  close  in  value  and  the  seasonal  and 


K variation  of  the  parameters  is  very  similar.  The  solar-cycle  dependence 
P 

of  foF2  appears  to  be  adequately  described  by  the  ITS  model  since  there  is 


very  little  spread  in  the  values  of  A for  a given  month.  The  width  and 

m 

center  position  along  with  other  data  were  used  to  determine  the  position  of 
the  southern  wall  of  the  auroral  oval  as  a function  of  corrected  geomagnetic 
time. 

We  have  combined  all  of  these  results  in  a set  of  corrections  to  be 
applied  to  the  median  model.  In  our  previous  work  based  on  the  1958  data, 


we  also  noted  an  apparent  linear  dependence  of  the  K correction  during 

P 


the  solar  day.  In  that  data,  there  also  seemed  to  be  a seasonal  dependence 


TABLE  1-2  High-Latitude  Trough  Parameters 
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of  this  correction.  The  data  for  the  other  years  does  not  show  this  effect. 

The  K correction  was  determined  to  be 
P 

A (K  ) = -0. 10  + 0. 025  N 
P gP 

ITS 

and  foF2A  = foF2  (1  - A (K  )) 

A p 

This  correction  of  between  -7. 5%  and  +7.  5%  is  applied  independently  of 
space  and  time  to  the  model. 

The  spatial  dependent  features  of  the  trough  and  auroral  oval  have 
been  taken  to  be  of  the  following  form: 


f°F2  = foF2A  (1  + An) 

where  AN  = Aq  Jexp  XA  exp(-XA2/2)  f°r  & 0A 

for  the  auroral  oval 


0A  = 71.9  - 2.5K  - 5.1  cos  I ^ (t  -1)"1  (degrees) 

A p L12  cgm  J 

0 s 0 s 0 

mn  A mx 

0 = 68.9  - K - 5.  1 cos  T "rr  (t  -1)1 

mn  p LI  2 cgm  J 

0 = 70.9  - K - 5.  1 cos  [" (t  -I)] 

mx  p L12  v cgm  J 
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. *4+ 


7.0 


(degrees) 


4°  s XT  s 6° 


X = (X  - 0A)/X 

g m AT 


XA  = 1.5  X 
A g 


and  for  the  trough 


* ■ T [I  * ».  HH] 


for  X < 0 
m A 


where 


T = Tj  7exp  X exp  (-  X^/2)  exp  (-(tcgm-3)2/12)  for  |Xg|sl 


T = Tj  exp  (-2.5  (X  -1)  ) for  |Xg|>l 


= 0 for  Xs90 


or  6. 0 < t <18.0 
cgm 


T X > 94. 6 and  t <6.0  or  t ^ 18.0 
o cgm  cgm 


U 

^ ~ 90°  $xs  94.  6 and 


° 4.6° 


t < 6.  0 or  t & 1 8.  0 
cgm  cgm 


and  where  t is  the  corrected  geomagnetic  time.  The  maximum 
cgm 

amplitudes  have  been  taken  to  be  T = -0.  20  and  A =0.  20.  Examining 
r o o 


■T  ’■■■'&  i > . . >.  . . 
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A 


the  five  years  data  lead  to  the  conclusion  that  only  a small  variation  of  the 

trough  position  and  width  should  be  contained  in  the  corrections.  This  is 

contained  in  the  above  corrections  by  setting  maximum  and  minimum 

values  for  0 and  X . The  limits  on  X and  the  functional  form  used 
A 1 ! 

allow  the  full  width  at  one  fourth  maximum  of  the  trough  region  to  be  be- 

o o 

tween  o.  3 and  9-  5 . 

The  resulting  model  for  foF2  was  used  to  calculate  the  RMS  residual 
for  the  same  grouping  of  data  as  was  done  with  the  ITS  model  that  was 
shown  in  Figure  I- 1 and  1-2.  Grouping  the  data  according  to  solar  day  and 
night  and  month,  24  values  of  were  calculated.  The  number  of  months 
in  which  the  new  model  showed  an  improvement  over  the  ITS  model  is 
shown  in  the  following  table. 


Year 

Solar  Day 

Solar  Night 

I960 

10 

1 

1962 

11 

10 

1964 

10 

8 

1966 

9 

10 

1968 

10 

10 

The  use  of  the  model  depends  upon  having  some  prior  knowledge  of  the 
activity  in  the  polar  region.  During  periods,  such  as  1958,  in  which  the 
depression  of  foF2  due  to  the  trough  is  particularly  weak  a smaller  value 


of  the  amplitude  should  be  used.  When  the  model  was  applied  to  a 

small  sample  of  polar  stations  measurements  of  foF2  made  during  19^9 

and  1970,  no  significant  improvement  over  the  ITS  model  was  made 

when  a T of  20%  was  used, 
o 
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II.  STUDIES  OF  IONOSPHERIC  RAY  TRACING  IN 
PARABOLIC  IONOSPHERES  WITH  GRADIENTS 


INTRODUCTION: 

A fundamental  problem  in  the  design,  planning,  and  operation  of  HF 
communications  systems  operating  in  the  frequency  range  2 to  30  MHz 
is  to  obtain  realistic  quantitative  estimates  of  the  expected  long  term 
performance  of  such  systems.  Such  predictions  are  usually  effected  by 
ray  tracing  techniques  through  a realistic  model  of  the  ionosphere.  In 
the  general  case,  where  realistic  global  electron  density,  magnetic  field, 
and  electron  collision  frequency  models  are  taken  into  account,  numerical 
integration  of  the  Haselgrove  equations,  ^ such  as  described  by  Jones^*, 
is  generally  required.  Such  an  approach  generally  requires  powerful 
computer  facilities  beyond  the  capability  of  field  operations. 

Other  techniques  generally  employ  a spherically  symmetric  electron 

density  model  neglecting  the  magnetic  field.  In  such  an  ionosphere  Bouguer's 

rule  applies  and  the  ray  trace  equations  may  be  reduced  to  quadratures  (see, 

(3) 

for  example  Croft  and  Hoogasian  ).  These  integrals  have  approximate 
solutions  for  the  parabolic  electron  density  model  which  were  used  by 

(4) 

Barghausen  et  al  to  develop  the  ITS- 78  operational  parameter  prediction 
program. 

For  certain  simple  ionospheres  these  integrals  have  exact  solutions; 
the  solution  for  the  quasiparabolic  ionosphere  is  described  by  Croft  and 
Hoogasian^  \ the  quasilinear  case  by  Westover 

These  solutions,  based  on  a spherically  symmetric  ionosphsre,  may 
be  employed  for  cases  where  horizontal  gradients  exist  in  electron  density 
if  the  effect  of  such  gradients  may  be  adequately  represented  by  a displace- 
ment of  the  center  of  curvature  of  the  spherically  symmetric  ionosphere 
from  the  center  of  the  earth.  The  formal  solution  in  such  a situation  is 
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described  by  Folkestaa  . A method  of  determining  the  effective  center 
of  curvature  based  on  variation  of  virtual  height  was  developed  by 

(7)  ^g j 

Beckwith'  as  a basis  of  predicting  propagation  characteristics' 

Fundamental  to  all  HF  computer  prediction  programs  is  a synoptic 

numerical  representation  of  the  ionospheric  characteristics.  Typically, 

electron  density  models  are  developed  by  representing  the  ionosphere  as 

one  or  more  parabolic  layers,  the  parameters  of  which  (critical  frequency, 

layer  height  and  semithickness)  are  determined  from  global  maps  of  co- 
(4) 

efficients  or  empirical  functions  representing  diurnal,  seasonal,  and 
solar  cycle  variation.  Implicit  in  these  models  are  horizontal  gradients 
in  electron  density  due  to  variations  in  the  layer  parameters.  It  is  useful 
to  inquire  into  the  effect  of  such  gradients  in  the  propagation  characteristics 
of  HF  radio  waves.  In  particular,  it  is  useful  to  inquire  as  to  whether  the 
effects  of  such  gradients  may  be  reduced  to  a displacement  of  the  center  of 
curvatire  of  an  equivalent  spherically  symmetric  ionosphere.  If  so,  the 
formalism  exists  to  develop  a rapid,  efficient  and  accurate  algorithm  for 
the  prediction  of  propagation  characteristics. 

The  present  study  is  an  investigation  of  the  effects  of  gradients  in  the 
critical  frequency.  Numerical  ray  traces  were  calculated  using  the 
ARCON  II  3D  ray  tracing  program;  a dipole  approximation  of  the  geomagnetic 
field  was  included  with  a simple  parabolic  electron  density  to  represent  the 
F2  layer.  Ad  hoc  gradients,  parallel  and  perpendicular  to  the  ray  direction, 
were  included.  The  results  of  the  calculation  were  compared  to  those  de- 
rived from  the  analytic  solution  of  the  quasiparabolic  electron  density  model. 
The  rotation  of  the  rays  was  investigated  in  its  dependence  on  propagation 
characteristics,  particularly  the  virtual  height. 


COORDINATE  SYSTEMS  AND  MAGNETIC  FIELD: 


The  ARCON  II  3D  Ray  Tracing  Program,  as  its  predecessor  the 
JONES  3D  Ray  Tracing  Program,  explicitly  considers  two  coordinate 
systems:  (1)  The  geographic  coordinate  system  with  the  pole  aligned 

with  the  axis  of  rotation  of  the  earth  and  prime  meridian  passing  through 
Greenwich,  England,  and  (2)  the  geomagnetic  coordinate  system  with  pole 
aligned  with  the  dipole  magnetic  pole,  and  prime  meridian  containing  the 
geographic  north  pole.  Input  and  output  coordinates  are  generally  refer- 
enced to  the  former,  while  internal  calculations  are  referenced  to  the 
latter. 

For  ray-tracing  applications  wherein  the  geomagnetic  field  is  ignored, 
the  differences  between  the  two  is  moot;  similarly,  when  a simple  isotropic 
model  is  employed  for  the  electron  density.  However,  when  ray-tracing 
with  magnetic  field  and  a realistic  global  electron  density  model,  ARCON  II 
constrains  the  user  to  reference  the  electron  density  model,  collision 
frequency  model  and  geomagnetic  field  model  to  a single  coordinate  system. 
Specifically,  the  geomagnetic  field  model  is  the  IGRF  (1965)^'’  model 
referenced  to  the  geomagnetic-dipole  coordinate  system;  thus,  the  electron 
density  must  be  referenced  to  this  system. 

It  is  often  desirable  for  the  computational  coordinate  system  to  be 
specified  arbitrarily;  in  particular,  ray  tracing  applications  ultimately 
require  range  and  azimuth  of  the  salient  features  from  the  transmitter. 
Range  and  azimuth  are  the  spherical  coordinates  of  a coordinate  system 

with  pole  at  the  transmitter  and  prime  meridian  containing  the  north  geo- 

* 

graphic  pole,  this  coordinate  system  shall  be  referred  to  as  the  trans- 
mitter coordinate  system  (TCS),  and  is  patently  the  natural  coordinate 

system  fcr  the  general  ray  tracing  problem. 

* 

In  a right-handed  coordinate  system,  the  "longitude"  coordinate  is 
II  - azimuth. 
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One  application  where  the  use  of  the  TCS  is  particularly  desirable 
is  in  the  present  study  of  the  effects  of  longitudinal  and  transverse 
gradients  of  the  electron  density  on  propagation  characteristics.  The 
TCS  is  the  canonical  coordinate  system  for  this  study;  the  desired  pertur- 
bations to  the  electron  density  model  may  be  modelled  directly  and  simply; 
the  resultant  effects  on  propagation  characteristics  are  also  simply 
presented. 

Another  example  is  in  using  the  numerical  electron  density  model, 
the  default  model  of  ARCON  II.  Generally,  a transmit  site  and  coverage 
domain  are  predefined  for  the  application,  and  a numerical  ionosphere  is 
generated  sufficiently  extensive  to  trace  the  required  rays.  Often,  the 
required  application  is  stated  as  a transmitter  at  a given  geographic 
location  covering  a small  azimuth  domain  to  several  thousand  kilometers 
in  range.  In  constructing  the  numerical  ionosphere,  the  user  must  convert 
this  information  to  an  equivalent  domain  in  geomagnetic  coordinates. 

Usually,  the  desired  coverage  requires  a considerable  larger  domain  in 
geomagnetic  coordinates  than  in  the  TCS;  a considerable  portion  of  the 
resulting  array  of  plasma  frequencies  remains  inactive  during  the  ray 
trace,  and  the  user  may  be  constrained  to  use  a larger  grid  spacing  than 
he  would  otherwise  deem  desirable  in  order  to  limit  core  storage  require- 
ments.  Use  of  the  TCS  in  the  ionosphere  construction  program  precludes 
the  requirement  on  the  user  to  transform  coordinates  and  allows  for  an 
ionosphere  tailored  to  the  application. 

In  conclusion,  we  submit  that  the  TCS  is  the  preferred  coordinate 
system  in  most  ray-tracing  applications. 

However,  when  the  coordinate  system  for  the  ray  tracing  algorithm 
is  arbitrarily  specified,  there  remains  the  problem  of  the  proper  treatment 
of  the  geomagnetic  field.  The  present  geomagnetic  field  model  is  the  IGRF 

| 
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(1965)  model  referenced  to  the  geomagnetic  dipole  coordinate  system.  It 
is  undesirable  in  terms  of  execution  efficiency  to  effect  a coordinate  trans- 
formation to  the  coordinate  system  of  the  geomagnetic  field  model  every 
time  the  field  components  are  desired.  The  solution  of  the  problem  is  to 
rotate  the  components  of  the  geomagnetic  field  and  reference  the  field 
directly  to  the  chosen  computational  coordinate  system. 

The  geomagnetic  field  vector  B is  given  by 

^ = - grad  V,  (1) 

where  V is  a scalar  potential  given  by 

B 

V(r,  0 , cp)  = S V (r,0,cp)  (2) 

„ n 
n=0 

The  partial  potentials  V are  given  by 


,r  . to  , , .n-t-1  rn.  ( m , ,m  . 1 

V (r,0,cp)  = r (r  /r)  P (cos0)4g  cosmep  + h sinmcpi 

n 0-0  n In  n j 


(3) 


m=0 


where,  r is  the  radius  of  the  earth,  r,0,cp  are  the  coordinates  of  the 
o 

field  point,  the  g and  h are  the  expansion  coefficients  given  by 
n n 

1AGA  commission  2 in  geographic  coordinates  and  equivalently  by 
Mead'1”'  in  geomagnetic  dipole  coordinates,  and  the  are  the  Schmidt 
quasi-normalized  associated  Legendre  functions  of  degree  n and  m 
given  by 


Pm(X)  = 

n 


e (n-m) ! 
m 

_ (n+m) ! 


1/2  m/2  m+n 

l.  d 

(1-X  ) 


dx™ 


+ n 


(x2  -i)n 


(4) 
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where  e = 2 for  msl. 
m 


Consider  the  dipole  terms  (n=l)  of  equation  (3): 

V^r,  0, CP)  = (rQ/r)3  (g J x + h*  y + z)  (5) 

where  x = r sin  0 cosep 
y = r sin  Q sincp 
z = r cos  0 

Equation  (5)  may  be  expressed  in  vector  notation  as 

Vj(?)  = (rQ/r)3P.?  (6) 

where  the  components  of  the  dipole  moment  vector  ^ are  identified  as 


<P*-'Py’P,>  ' 


Given  the  components  of  the  dipole  moment  ^ in  the  coordinate 
system  r,  then  the  components  in  a rotated  coordinate  system 


where  3 is  the  rotation  matrix  transforming  between  the  two  coordinate 
systems.  Equation  (9),  then,  is  the  desired  transformation  of  the  dipole 
field  to  an  arbitrary  coordinate  system. 

The  higher-order  components  of  the  geomagnetic  field  are  elements 
of  higher- rank  tensors  and  may  be  transformed,  in  principle,  by  appli- 
cation of  the  rotation  matrix.  However,  it  may  be  more  efficient  compu- 
tationally to  take  advantage  of  the  scalar  nature  of  the  and  generate 

2n+  1 simultaneous  equations  for  the  n^  order  expansion  coefficients. 

(10) 

The  latter  approach  was  used  by  Mead  to  develop  the  values  of  the 
components  in  the  dipole  coordinate  system. 

In  the  ray-trace  applications  currently  studied,  the  dipole  approx- 
imation of  the  geomagnetic  field  is  adequate;  thus,  only  the  transfor- 
mation of  the  dipole  terms  given  by  equation  (9)  has  been  implemented 
for  use  in  the  ARCON  II  ray  trace  program.  However,  we  hope  to 
optimize  the  versatility  of  ARCON  II  in  the  near  future  by  implementing 
an  algorithm  to  rotate  the  higher  order  components  of  the  geomagnetic 
field. 


ELECTRON  DENSITY  MODEL: 

The  ambient  ionosphere  for  this  study  was  specified  as  a parabolic 
electron  density  model.  The  plasma  frequency  f^  in  this  model  is  given 


f*  = *1  (i-y2)  h -y<h<h  +y 
pc  m m 


f = 0 
P 


othe  rwise 


where 


y = (h-h  )/y 

m m 


and  f^  is  the  critical  frequency,  the  parabolic  layer 

height  of  maximum  electron  density,  y^  the  layer  semithickness,  and 

h the  altitude.  For  this  study,  h and  y were  fixed  at  300  km  and 

m m 

100  km  respectively,  and  the  critical  frequency  f was  taken  nominally 

to  be  10  Mhz.  However,  ad  hoc  gradients  were  allowed  in  f such  that 

c 


fc  = fco  + A 0O)  + B sin  e (cp-c pq)> 


where  f = 10  MHz,  A and  B are  ad  hoc  gradients,  0 is  range  from 

the  transmitter  and  Cp  is  the  azimuth  of  the  ray  point  from  the  receiver. 

Since  the  actual  variability  of  fQF^  is  generally  less  than  10  MHz 

over  30  degrees,  values  of  A and/or  B less  than  2 (with  f =10  MHz 

co 

and  0,  cp  in  radians)  suffice  to  simulate  typical  ionospheric  conditions. 
Sets. of. ray  traces  were  computed  for  A and/or  B equal  to  0,  1,  and  2. 


I 


Since  it  was  desired  that  the  nominal  electron  density  profile  (i.  e. 

f =10  MHz.)  obtain  at  or  near  the  ray  reflection  point,  the  reflection 
c 

point  (9^,  cpQ)  must  be  estimated  prior  to  tracing  each  ray.  Since  the 
azimuthal  deviations  were  small  (generally  less  than  2 degrees),  cp^ 
may  be  taken  as  the  initial  takeoff  azimuth.  The  nominal  range  0^  was 
then  estimated  from  virtual  height  analysis  appropriate  to  parabolic 
layer  theory. 

The  ray  geometry  is  shown  in  Figure  3.  Parabolic  layer  theory 
gives  the  virtual  height  h^  at  operating  frequency  f and  takeoff  elevation 
as. 


wh  e re , 


and, 


, , 1 . , 1 + u, 

h = h + — y u-tog  (- 

v b 2m  1-u 


h = h - y 
b m m 


u = (f/f  ) sinB 
c v 


(1?) 


(13) 


while  3^  is  the  local  elevation  angle  of  the  unrefracted  ray  at  the  apex 
of  the  virtual  triangle  and  is  given  by. 


cosg  = (r  /r  ) cosB 
Kv  o v o 


(14) 


where  r is  the  radius  of  the  earth,  and  r = r + h . 

o v o v 

Equations  (12)  to  (14)  patently  require  iteration  for  the  solution 

h ; starting  from  the  initial  estimate  h - h,  two  iterations  were 
v v b 


considered  sufficient  to  estimate  the  range  0q  to  the  reflection  profile 
given  by, 

0 = B - P 

o pv  o 

This  technique  located  the  reflection  profile  to  better  than  4 
kilometers  in  most  cases  for  the  control  ionosphere  (no  gradients 
and  no  magnetic  field).  Larger  errors  only  occurred  at  elevations 
close  to  the  skip  distance  elevation;  however,  all  deviations  were 
less  than  16  kilometers  at  elevations  less  than  the  skip  elevation. 

When  gradients  and  magnetic  field  were  introduced,  the  actual 
points  at  which  reflection  occurred  never  deviated  from  the  estimated 
profile  by  more  than  100  kilometers;  the  deviations  exceeded  32  kilo- 
meters only  near  the  skip  distance  elevation  at  frequencies  exceeding 
25  MHz.  Thus,  the  critical  frequency  at  the  reflection  point  never 
differed  from  the  nominal  value  of  10  MHz  by  more  than  5%,  and 
generally  by  much  less  than  1%. 
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THE  QUASIPARABOLIC  ELECTRON  DENSITY  MODEL  (QP  MODEL) 


The  QP  Model  has  considerable  utility  as  an  initial  approximation  to 
the  ray  traces  in  the  ionospheres  under  consideration,  primarily  because 
it  yields  analytic  expressions  for  ray  traces  through  this  model  ionosphere 
as  shown  by  Croft  and  Hoogasian.  ^ The  QP  Model  has  been  used  by  Rao*111 
to  develop  a method  of  extracting  ionospheric  parameters  from  measured 
minimum  group  paths  at  a set  of  frequencies. 

The  QP  Model  gives  the  plasma  frequency  f^( r)  in  terms  of  ionospheric 

parameters  of  a single  layer  ionosphere,  to  wit,  critical  frequency  f^, 

layer  altitude  h , and  semithickness  y . 

’ m 7m 


f (r)  = f fl  - y2~l  2 , r < r < r r /(r  - y ) 

p'  ' cL  J b mbb  'm 


f (r)  = 0 otherwise 
P 


(16) 


where  y = (r  / r).(r-r  )/y  ’ rm  = r0  + V 


r,  = r -y  , and  r„  is  the  radius  of  the  earth, 

b m m 0 


The  ray  trace  solutions  derived  by  Croft  and  Hoogasian  are  shown 

as  families  of  group  path  and  range  as  functions  of  takeoff  elevation  in 

Figures  1 and  2 respectively.  The  parameters  of  the  QP  Model  used 

were  r = 6370  km,  f = 10  MHz,  h = 300  km,  and  y = 100  km. 
0cm  m 

Each  curve  is  for  a constant  f,^/fc  where  the  transmitter  frequency, 

f^,,  varies  between  2 and  30  MHz  in  2 MHz  steps. 
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Our  objective  is  to  use  knowledge  of  the  ionosphere  in  lieu  of  precise 
knowledge  of  elevation  angle  to  determine  ground  range  from  a measured 
group  path  at  a given  frequency.  Thus,  we  need  to  enter  Figure  1 at  a 
given  group  path  and  frequency,  extract  an  elevation  angle  and  enter 
Figure  2 for  range;  in  short,  we  desire  to  invert  the  functions  of  Figure  1. 

To  do  so  efficiently,  we  analyzed  the  solutions  of  Croft  and  Hoogasian, 
simplified  the  functional  forms,  and  chose  a convenient  variable  with 
which  to  represent  the  solution. 

The  ground-to-ground,  one-hop  group  path  P'  is  given  by 

P'  = 2r  sinY  (G  + G ) (17) 

b p r 1 


where  r^  (previously  defined)  is  the  base  of  the  ionosphere,  Y^  is  the 
maximum  entry  angle  above  which  all  rays  at  that  frequency  will 
penetrate  the  ionosphere,  G is  the  contribution  from  the  free  space 
propagation  (transmitter  to  bottom  of  ionosphere),  and  G^  is  the 
contritution  from  integration  within  the  ionosphere  (bottom  of  iono- 
sphere to  ray  apogee).  The  contributions  G and  G are  giver,  by: 

JC  1 


I 


The  constants  G^,  C ^ and  siny^  depend  upon  the  transmitter 
frequency  f and  are  given  by 

Gj  = F2/A 
G2  = B/2B' 

siny  = B' IT  Ja  (20) 

P 

where, 


Note  that  at  F 


= F , v =0;  thus,  F is  the  maximum  (scaled) 
m 'p  m 

frequency  for  which  the  QP  model  ionosphere  will  support  reflections. 
However,  in  ground  to  ground  transmission  y^  can  be  no  less  than  y^; 
accordingly,  the  maximum  (scaled)  reflecting  frequency  ^max  *s  con- 
siderably less  than  F , and  is  determined  from  y.,  which  yields, 

1 m 0 


F 

max 


(2-) 


where 


a = 1 + 7 (F  tany  )2 
c.  m U 

b = secy^ 


(^ 


To  determine  the  minimum  group  path  at  a given  frequency,  it  is 
necess  ary  to  solve  the  equation. 


G =dG/3z 
z 


1 - z/z' 


2G1G2/(1-z  ) -Gl 


= 0 


(24) 


Although  equation  (24)  reduces  to  a cubic  in  z , thus  permitting  a 
solution  in  closed  form,  it  is  as  effective  numerically  to  employ  a 
Newton  iterative  method,  for  which  the  second  derivative  of  G is 
required,  and  given  by. 


G = d^G/dz2  = ((z/z1)2  - 1 ) / z 1 + 4G  G z/(l  - z2)2 

ZZ  1 L 


(25) 
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Solutions  of  equation  (24)  obtain  for  all  F<Fmax  at  z between  zQ  and 
unity,  favoring  the  upper  limit.  Thus,  the  iterative  procedure  is 
started  from 


0.8  + 


0.  2 z 


0 


and  the  solution  z taken  as  the  limit  of  the  series 
m 


(26) 


z(n)  = z(n_1)  - G /G  . (27) 

mm  z zz 


With  the  salient  features  of  the  group  path- elevation  curve  at  a 
given  frequency  thus  established,  equation  (17)  may  now  be  inverted 
for  elevation,  given  group  path.  A unique  solution  exists  in  z^  <z<z^ 
if  the  input  frequency  is  less  than  F , and  the  given  group  path  is 
between  the  minimum  group  path  and  the  group  pp.ch  at  z = z().  Since 

equations  (18)  are  intrinsically  irrational,  an  iterative  technique  is 
mandated.  From  Figure  1,  it  is  clear  that  a good  first  estimate  z^ 
to  the  low  angle  solution  can  be  obtained  by  first  subtracting  from  the 
given  group  path  the  excess  AP  over  specular  reflection  from  the  base 
level  at  elevation  g = 0,  and  treating  the  remainder  as  specular 
reflection  (see  Figure  4.).  Therefore,  if  P'  is  the  given  group  path  at 
frequency  f, 

then  write  P1  s 2rQsinYp  ^g  + Gj(zo)j 

G = P'/2r  siny 
p op 


i 

i 


3 
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Fig.  4.  Initial  Estimation  for  Iterative  Inversion  of  Group 
Path  Length  vs.  Elevation  Curve  of  Figure  1. 
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R * Gf(z(0)) 


from  which  z^  ^ = T ^ Zo 1 


If  the  high-angle  solution  is  desired,  the  initial  estimate  may  be 
taken  between  that  for  minimum  group  path  and  penetration: 


z(0)  = (1  + z )/2 

m 


The  solution  Z is  taken  as  the  limit  of  the  series 


z(n-I)trc_ 


J z 


From  this  solution  the  entry  angle  V,  takeoff  elevation  |3,  ground 
range  D and  penetration  depth  L are  given  by 


. -1, 

sin  (z  siny  ) 
P 


cos  (cosy/cos  Y 


= 2r  f (Y-P)  + (cos  y/Jac)  -tog 

O L 

= rb(B'/A)  [l  - yl-zj 


where 


C = B + B'  + (FzsinY  ) 

P 


= [l  + TcTa  z + (B'/A)  z*]  /T 1 


[ 1 - Z'] 


DATA  FILES  GENERATED  BY  THE  NUMERICAL  RAY  TRACE 


Six  files  containing  results  from  the  ARCON  II  3D  ray  tracing  program 
for  frequencies  from  2 to  30  MIIz  and  elevations  2*  to  45°  (or  penetration) 
were  c reated  for  further  statistical  study.  Parameters  kept  fixed  for  all 
files  are  listed  in  Table  I,  while  Table  II  lists  elements  that  were  changed 
from  file  to  file.  A record  was  generated  for  each  ray  traced  consisting  of 
the  basic  coordinate  data  of  the  ray  at  the  following  points:  transmitter, 

entry  to  ionosphere,  ray  apogee,  exit  from  ionosphere,  and  ground  re- 
flection. In  addition,  information  identifying  the  ray  such  as  frequency, 
elevation,  and  information  identifying  the  ionosphere  was  appended. 


Table  1:  Fixed  Ray  Trace  Parameters,  All  Files 


TRANSMITTER  LATITUDE: 
TRANSMITTER  LONGITUDE: 
TAKEOFF  AZIMUTH: 
ELEVATION  RANGE: 
ELEVATION  STEP: 
FREQUENCY  RANGE: 
FREQUENCY  STEP: 
INTEGRATION  STEP  SIZE: 


45°  N 
68?  W. 

045° 

Cf  to  45°,  or  penetration 
1° 

2 MHz  to  30  MHz 
2 MHz 
0.  25  km 
Parabolic 


MODEL  IONOSPHERE: 

NOMINAL  CRITICAL  FREQUENCY:  10  MHz 

300  km 


LAYER  HEIGHT: 
SEN  ITHICKNESS: 


100  km 


TABLE  II:  VARIED  RAY  TRACE  PARAMETERS 


FILE 

FIELD 

LONGITUDINAL 

TRANSVERSE 

NUMBER 

MODEL 

GRADIENT  (A) 

GRADIENT  (B 

1 

NONE 

0 

0 

2 

DIPOLE 

0 

0 

3 

DIPOLE 

1 

0 

4 

DIPOLE 

2 

0 

5 

DIPOLE 

0 

1 

6 

DIPOLE 

0 

2 

1 


ANALYSIS  OF  RESULTS: 


The  results  of  ray  traces  through  the  ionospheres  with  longitudinal 
gradients  in  the  critical  frequency  were  analyzed  to  discover  correlations 
between  deviations  in  range  and  group  path  from  the  predictions  of  the 
unperturbed  quasiparabolic  electron  density  model. 

The  control  set  of  ray  traces,  no  field  and  no  gradients  in  critical 
frequency,  yielded  results  from  the  numerical  solution  of  the  Haselgrove 
equations  substantially  in  agreement  with  the  analytical  solution  of  the 
quasiparabolic  model;  the  root  mean  square  deviations  of  the  differences 
in  range  and/or  group  path  length  between  the  two  were  of  the  order  of  the 
integration  step  size  (0.  25  kilometers). 

Tie  essential  effect  of  turning  on  the  (dipole)  magnetic  field  was  to 
introduce  a certain  amount  of  "noise"  in  the  ray  trace  results  although 
systematic  frequency  dispersion  is  apparent.  At  frequencies  exceeding 
two  megahertz,  the  root  mean  square  deviations  in  range  and  group  path 
from  the  QP  solution  were  less  than  11  kilometers.  At  two  megahertz, 
the  RMS  deviation  was  less  than  20  km.  Figure  5 shows  the  range  deviations 
as  functions  of  group  path  length.  The  subsets  for  frequencies  2,  4,  6 MHz, 
etc.  are  indicated  by  characters  B,  D,  F etc.  Significant  range  deviations 
occur  primarily  at  low  frequencies  (2  and  4 MHz,  when  the  gyrofrequency  is 
significant  compared  to  the  operating  frequency),  and  at  other  frequencies 
near  the  minimum  group-path  length. 

The  additional  effects  of  longitudinal  gradients  in  the  critical  frequency 
are  shown  in  Figure  6 (A  = 1)  and  Figure  7 (A  = 2).  No  significant  effect 
is  noticed  in  either  case  at  2 End  4 megahertz;  at  higher  frequencies  the 
gradient  effect  tends  to  dominate  over  the  magnetic  field  effect.  A definite 
frequency  dispersion  is  evident  in  these  range  deviations  although  patently 
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Range  Deviations  from  Quasiparabolic  Solutions  vs.  Quasiparabolic 
Group  Path  Length.  Case  (a).  Dipole  Field-no  Gradients. 
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Fig.  6.  Range  Deviations,  Dipole  Field-longitudinal  Gradient  (A-l) 


Group  Path  (km) 


Group  Path  (km) 


nonlinear  in  total  group  path  length.  The  range  and  group  path  deviations 
from  the  quasiparabolic  model  were  analyzed  in  an  attempt  to  construct  a 
model  of  the  following  form 


R = RQp  + a (f.  A)  + b (f.  A)  (P-  - P'Qp) 


where  R and  P'  are  ground  range  and  group  path  length  as  calculated  from 
the  numerical  ray  trace,  RqP  and  P'qp  are  calculated  from  the  quasi- 
parabolic model,  while  a (f.  A)  and  b (b.  A)  are  functions  of  frequency 
and  longitudinal  gradient  of  critical  frequency  f 


1 Sfc 

fc  a© 


Table  III  lists  the  results  of  this  statistical  analysis  giving  th<  fitting 
parameters  a and  b together  with  the  standard  deviations  a for  A = 1 and 
2.  A moderately  well  defined  frequency  dispersion  in  a (f)  is  readily 
apparent  (Figure  8)  and  appears  to  be  directly  proportional  to  A.  On  the 
other  hand  b (f)  appears  not  to  be  a well  defined  function  of  frequency,  but 
statistically  independent  of  A. 

Figures  9 and  10  show  the  range  deviations  as  functions  of  group  path 
length  for  transverse  gradients  in  critical  frequency  (B  = 1 and  2 respectively), 
while  Figures  11  and  12  present  the  corresponding  azimuthal  deviations. 

Each  pair  reinforces  the  conjecture  that  gradient  effects  are  directly  pro- 
portional to  the  magnitude  of  the  gradient.  In  addition,  range  deviations 
appear  to  be  strongly  correlated  with  azimuthal  deviations  while  both  are 


M * 
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TABLE  III:  Fit  of  range  deviations  to  group  path 
deviations  for  the  ionospheres  with 
longitudinal  tilt. 


A=  1 


A=2 


A=0 


Fig.  8.  Frequency  Dispersion  in  Range  Deviation  Fitting  Parameter. 


CD  Azimuthal  Deviation 
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Azimuthal  Deviation 


extremely  sensitive  to  small  changes  in  elevations  near  the  minimum 
group  path  length.  Furthermore,  frequency  dispersion  of  the  transverse 
effects  is  well  defined. 

Figures  13  and  14  present  azimuthal  deviations  as  functions  of  takeoff 
elevation. 

Several  interesting  features  of  the  effects  of  gradients  in  critical 
frequency  are  worthy  of  notice.  In  the  first  place,  effects  due  to  the 
magnetic  field  are  small  and  important  primarily  at  lower  frequencies. 

This,  of  course,  was  expected  since  the  gyrofrequency  is  of  the  order  of 
1 megahertz.  Secondly,  the  effects  exhibit  a well  separated  frequency 
dispersion.  Thirdly,  the  effects  appear  to  be  fundamentally  proportional 
to  the  magnitude  of  the  appropriate  derivative  of  critical  frequency. 

It  does  not  appear  likely  that  these  features  are  to  be  simply  related 
to  parameters  of  the  total  path  such  as  total  group  path  length  and  ground 
to  ground  range.  These  quantities  include  considerable  free  space  propa- 
gation, more  or  less  depending  on  takeoff  elevation.  The  important  contri- 
bution for  this  study  is,  of  course,  from  the  integration  through  the  iono- 
sphere. 

A productive  method  of  fast  and  accurate  ray  tracing  is  to  use  a simple 
electror  density  model,  such  as  the  quasiparabolic  model,  which  has  an 
analytic  solution  to  the  ray  trace  equations,  and  allow  this  ionosphere  to  be 
tilted;  i.  e. , the  center  of  curvature  of  the  spherically  symmetric  ionosphere 
is  displaced  from  the  center  of  the  earth.  Such  a "tilted"  ionosphere  is  then 
described  by  the  layer  parameters,  and  additionally  by  parameters  specifying 
the  center  of  curvature  and  lift  of  the  model  ionosphere.  Clearly,  the  latter 
quantities  will  depend  on  the  gradients  of  the  ionospheric  parameters,  and 
also  on  the  frequency  and  elevation  of  each  ray. 
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Azimuthal  Deviation  (degrees) 
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Fig.  13.  Azimuthal  Deviations  vs.  Elevation  (B=2). 
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Azimuthal  Deviation 
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Fig.  14.  Azimuthal  Deviations  vs.  Elevation  (B=l). 
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Consider  the  sets  of  ray  traces  through  the  ionospheres  with  longitudinal 
gradients  in  critical  frequency  (A  = 1 and  2).  Because  of  this  gradient,  a 
ray  entering  the  ionosphere  at  a local  elevation  angle  y will  exit  the  iono- 
sphere at  a local  elevation  angle  y'  4 Y*  The  effective  tilt  of  the  ionosphere 
for  this  ray,  then,  is  given  by  (y1  - y)/2. 

Inspection  of  a small  sample  of  these  ray  traces  revealed  that  the 
effective  tilt  is  strongly  correlated  with  the  penetration  depth  Ap  of  the  ray 
into  the  ionosphere.  This  correlation  for  all  rays  is  displayed  in  Figure  15. 
The  x axis  is  the  penetration  depth,  the  y axis  is  the  "normalized  effective 
tilt,"  NET  given  by 


NET  = (y'-y)/2A 


Plotted  points  "1"  are  rays  from  the  set  A = 1,  "2"  are  rays  from  the 
set  A = 2,  and  are  points  resulting  from  a least  square  fit  of  the  NET'S 

given  by 


NET  = a 


where  x = h„ly  > a = 1.  04  and  b = 1.37;  the  standard  deviation  of  the  fit 
P m 

for  all  rays  with  penetration  depth  exceeding  10  kilometers  (x  = 0.  1)  is  0.  11. 
This  fit  is  good  for  small  x but  fails  to  adequately  model  the  high  angle  rays. 
However,  the  most  important  feature  of  this  correlation  is  that  the  frequency 
dependence  is  nearly  completely  contained  in  the  penetration  death. 
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Virtual  Height 


A significant  improvement  in  analyzing  NET  can  be  effected  by  choosing 
the  virtual  height  h of  the  equivalent  ray  (see  Figure  3)  in  lieu  of  penetration 
depth  as  the  independent  variable.  The  resulting  correlation  is  shown  in 
Figure  16.  The  plotted  points  appear  to  be  levelling  off  for  the  high  angle 
rays;  three  terms  may  be  desired  to  improve  the  fit.  The  fitting  function 
plotted  is  given  by 


* 3 

NET  = a x + b x 


where  x = (h^  - h^/y^.  Other  functions  considered  to  model  NET  were 
a x,  and  a x + b x^.  The  fitting  parameters  and  standard  deviations  a of 
these  three  fitting  functions  are  listed  in  Table  IV. 


TABLE  IV:  Fitting  Parameters  for  NET 


Function 

a 

b 

0 

NET1  = ax 

. 606 

- 

. 039 

2 , 2 
NET  = a x + b x 

. 577 

. 045 

. 027 

* , 3 

NET  = ax  + bx 

. 598 

. 010 

. 033 

CONCLUSIONS: 


We  have  found  that  the  effective  tilt  of  a parabolic  ionosphere  with 
longitudinal  gradients  allowed  in  the  critical  frequency  can  be  well 
represented  by  a function  of  the  form 


T (f,  y 


af 

c 

3& 


(a  x + b x 


where 


T = (y'-Y)/2 

x = (h  -h  )/y 

v b m 

a = 0. 577  degrees  = 0.  01007  radians 
b = 0.045  degrees  = 0.00079  radians 


Angles  y1  and  y are  the  local  elevations  of  the  ray  at  exit  and  entry 
respectively  to  the  ionosphere  at  the  base  level  altitude  h^,  Yv  is  the 
elevation  angle  of  the  virtual  ray  at  the  virtual  height  h , f is  the 
operating  frequency,  f is  the  critical  frequency  of  the  layer,  and  y^ 
the  semithickness. 

The  displacement  D of  the  center  of  curvature  of  the  equivalent 
ionosphere  is  given  by 


D = (r  + h ) tan  T (r  + h ) T 

O V O V 


(33a) 


(33b) 
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This  result  leads  us  to  advance  the  following  conjecture  as  a basis  for 
formulating  a fast  ray  tracing  algorithm  in  ionospheres  with  gradients  in 
the  critical  frequency.  The  presence  of  gradients  in  the  critical  frequency 
of  the  reflecting  ionospheric  layer  may  be  accounted  for  by  assuming  an 
effective  tilted  spherical  ionosphere  the  center  of  curvature  of  which  is 
displaced  from  the  center  of  the  earth  by  the  vector 


D = -r  F (x)  grad  f /f 
f v c c 


(34) 


where 


F (x)  0.01007  x + 0.00079  xZ 


(35) 


and  the  gradient  is  to  be  evaluated  at  the  virtual  reflection  point. 

If,  on  the  other  hand,  the  critical  frequency  and  semithickness  did 
not  vary  but  gradients  existed  in  layer  height  then  the  ionosphere  is 
patently  tilted  with  the  center  of  curvature  given  by 


D = r grad  h (36) 

h v m ' ' 

Furthermore,  gradients  in  semithickness  alone  could  presumably 
be  reduced  to  an  equivalent  displacement  of  the  center  of  curvature 


D = r Y (x)  grad  y 
y v 6 7m 


(37) 
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Ln  addition,  the  equivalent  ionosphere  could  be  lifted  under  the 
influence  of  gradients  in  the  parameters  either  separately  or  because 
of  coupling.  This  "lift"  may  be  expressed  as 


= r^  L (x,  {v})  (38) 

where  the  argument  set  { V } includes  the  ionospheric  parameters  and 
their  gradients. 

No  conjectures  are  being  advanced  at  present  for  the  form  of  the 
functions  Y and  L. 

A possibly  useful  conjecture  to  consider  for  the  general  case  is 
that  the  center  of  curvature  of  the  equivalent  ionosphere  is  given  by 

D = Dr  + D,  + D + D.  (39) 

f h y v ' 
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III.  VERTICAL  ELECTRON  DENSITY  MODELS 
IN  THE  POLAR  REGION 


T 


In  Section  I,  we  have  described  improvements  made  to  the  ITS 
monthly  median  model  for  foF2  which  allow  for  a more  accurate 
desciiption  of  the  morphological  features  of  the  polar  ionosphere. 

In  earlier  work,  ^ we  have  described  models  applicable  in  the  polar 
region  for  some  other  ionospheric  parameters  such  as,  foE,  foFl, 
h E,  h E,  h FI  and  h F2.  These  parameters  were  combined  with 
several  empirical  models  to  produce  a three-dimensional  electron 
density  distribution  which  could  be  used  in  various  applications.  In 
one  application,  electron  density  distributions  are  obtained  on  a uniform 
spherical  grid  for  use  with  the  ARCON  Version  II  three  dimensional  ray 
tracing  program.  A grid  is  chosen  such  that  the  transmitter  is  at  the 
north  pole  of  the  coordinate  system.  This  coordinate  system  minimizes 
the  amount  of  data  needed  to  carry  out  any  ray  tracing  study. 

Extensive  studies  with  this  model  have  shown  some  areas  in  which 
the  ionosphere  is  not  correctly  represented.  We  will  describe  here  a 
new  formulation  for  a vertical  electron  density  distribution  when  one 
has  estimates  for  the  various  ionospheric  parameters  which  describe 
some  of  the  individual  layers.  One  area  in  which  an  improvement  was 
sough  was  in  the  E region  when  increased  ionospheric  activity  resulted 
in  an  elevated  value  of  foE  such  that  the  probability  of  there  not  being  a 
monotonic  vertical  electron  density  distribution  was  large. 

Since  the  ionospheric  parameters  for  the  various  layers  have  been 
determined  independently  of  each  other,  the  vertical  electron  density 
distributions  must  be  constructed  in  a manner  which  is  internally  con- 
sistent and  which  does  not  contain  any  spurious  horizontal  gradients. 
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r 


9 

An  example  of  an  inconsistency  would  be  a prediction  of  the  height 
for  an  FI  layer  at  a higher  altitude  than  that  for  the  maximum  density 
of  the  F2  layer.  In  general,  the  parameters  are  functions  of  the  geo- 
graphic and  geomagnetic  spatial  coordinates,  universal  time,  sunspot 
number  and  magnetic  activity.  The  following  Table  lists  the  parameters 
used  in  the  generation  of  a profile. 


Pa  rameter 

Needed  for  Profile 

Model  Used 

foF2 

R equired 

ITS  plus  modifications  (Section  I) 

h F2 
m 

R equired 

Reference  1 

y F2 
m 

R equired 

Reference  1 

foFl 

Optional 

Reference  1 

h FI 
m 

Optional 

Reference  1 

foE  (Solar) 
s 

Optional 

Reference  1 

foE  (Au  roral) 
A 

, h^E^  Optional 

Reference  1 with  modifications 

h E , y E 
ms  m 

Optional 

Reference  1 

The  parameters  for  the  FI  and  E region  are  in  general  functions  of 
local  time  (solar  zenith  angle)  and  may,  therefore,  be  present  only  in 
certain  geographical  regions.  We  have  attempted  to  allow  for  a smooth 
variation  of  the  profile  as  one  or  more  of  these  layers  appears  or  dis- 
appears from  the  ionosphere.  For  a given  set  of  conditions,  all  the 
desired  parameters  are  first  obtained  separately.  When  both  the  solar  E 
and  auroral  E layers  are  present,  these  parameters  are  combined  to 
produce  parameters  that  describe  a single  E layer.  These  two  para- 
meters, foE  and  h E are  generated  according  to  the  rules, 
m p 
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'Nf  * 


\ 


4 4 3 4 

(foE)  = (foEs)  + - (foEA) 


h E = (foE  )2  h E + (foE  J2h  E 
m sms  A m A 


(foE^)2  + (foE^)2 


The  electron  density  N(h)  above  the  maximum  height  of  the  F2 

layer  is  given  in  terms  of  two  profiles.  Between  h F2  and  a height  Z 

m c 

N(h)  is  defined  by  a Chapman  layer. 


N(h)  = F2  exp  [l  - § - exp  (- §)] 


where  § is  a scale  height, 


§ s h - h F2 
m 

y F2 
m 


N F2  = 1. 24  X 104  (foF2)2  el/cc. 
m 


The  altitude  Z - Z + h F2  the  low  end  of  the  topside  profile  is 
c t m r ^ 

obtained  as  the  solution  of  the  cubic  equation: 


Z + _CP_  Z + Y Z - y„  =0  where 
4 4 P 


the  coefficients  are 


and 


where 


F or 


where 

layer. 


P = -T 

Y = P2 


B 


1.  25  y F2 
m 


a = 3.14X10 


-4 


C = 5. 88  - 0. 15  X cos  X 

m 

r - 


1 2TT  , 

+ 0.41  cos  ( yj-r  (D  + 11);  sin  ^ - 0.14 


X = corrected  geomagnetic  latitude 
m 

\ = geographic  latitude 
X = solar  zenith  angle 
D = Day  number  (1-365) 

h s the  topside  profile  is 


N(h) 


N exp 
o r 


■^-C^tan  *(aCh) 


-1 

tan 


N = N(h  - Z + h F2)  as  determined  from  the  Chapman 
o t m r 
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Below  h F2,  there  are  four  possible  profiles  depending  upon  which 
layers  are  present  and  the  values  of  the  layer  parameters  relative  to 
each  other. 

If  no  E layer  is  present,  it  is  also  assumed  that  no  FI  layer  exists. 

In  this  case,  we  use  as  our  profile  model  a single  cubic  function  for  the 

electron  density  using  the  three  parameters  foF2,  h F2  and  y F2.  The 

m m 

thickness  of  the  layer  is  defined  as 


y 


A 


<Jl  y F2 
m 


The  electron  density  as  then  given  by 


N(h)  = N F2  I 1 + 2 [h  F2-h\3  - 3 (h  F2-h^2 
m / I m I I m 


Using  this  profile,  the  height  h = h F2  - 1.  3 y F2  is  such  that 
r ° m m 'm 


N (h  ) = 0. 01  N F2 
m m 


We  require  that  the  model  not  produce  a density  greater  than  this  value 

below  a specified  height  h . (taken  at  this  time  to  be  140  km. ) If 

mm 

h < h . , we  redefine  the  thickness  of  the  layer  such  that  N(h  ) = 
m min  min 

0.  01  N F2  by  the  following  shrinkage  of  the  layer: 


y 


A 


1.3 


(h  F 2 
m 


h . ) 

min 
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The  following  is  a sketch  of  a sample  profile  showing  the  relevant 
parameters . 


In  the  event  that  an  E layer  exists  but  that  there  is  no  FI  layer,  the 
profile  falls  into  two  categories.  In  the  first  case,  the  electron  density 
is  a monotonic  function  of  altitude.  In  the  second  case,  the  electron 
density  decreases  to  a minimum  value  between  the  heights  of  the  E layer 
and  F 2 layer  maximum  densities.  The  choice  of  which  profile  is  used 
is  based  upon  the  following  criteria. 

For  a monotonic  profile  to  be  used,  we  require  that  the  electron 
density  distribution  in  the  F2  region  not  fall  too  rapidly  so  as  to  be  in- 
consistent with  the  electron  density  at  the  E layer  maximum.  One 
possible  statement  of  this  condition  is  that 

h E ^ h F2  - Jz  y F2  (A) 

m m m 
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If  this  test  is  satisfied,  we  then  calculate  the  frequency  f^,  at  h^E  from 
a cubic  F2  layer  profile  using  the  following  expression: 


f 


2 

E 


(foF2)  1 + 2 


h F2  - h El 
m m 


The  second  condition  for  a monotonic  profile  is  that  f as  predicted 
from  a F2  layer  by  greater  than  f^E  the  value  determined  from  an 
independent  model  for  the  E region. 


f„  (h  E)  a f E (h  E) 
E m o m 


(B) 


When  this  condition  is  not  met,  the  F2  region  electron  density  has  dropped 
too  rapidly  to  be  consistent  with  an  independent  by  predicted  density  for  the 
E region. 

The  monotonic  profile  used  is  of  the  following  form: 


(i)  F2  region  h E^h^h  F2 

mm 


N(h) 


N F2  + (N  F2-N  E) 
m mm 


\h  I 
\ m 


-h  'l 

3 

^ h F2  - h 

m 

-h  E j 

i h F2-h  E 

m 1 

\ m m 
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(iii>  N(h)  = 0 


h < h E - ijz  y E 
m m 


This  profile  is  also  used  in  the  event  that  foE  ^ foF2. 

When  either  condition  A)  or  B)  is  not  satisfied,  a profile  of  the  type 
shown  in  the  following  figure  is  used. 


A parabolic  layer  is  used  in  the  F2  region.  This  profile  is  used  to 

determine  the  height,  h at  which  N(h  ) in  the  F2  region  is  equal  to 

r ii«  r L 

N E.  A cubic  equation  is  then  used  in  the  region  between  this  height 
m 

and  h E which  maintains  both  a continuous  electron  density  and  its 
m 

height  derivative.  If  we  take 


1 


then  we  require  the  following  conditions  for  the  coefficients 


(a) 

N(o)  = N E 
m 

= d 

(b) 

BN  , . 

— (°)  = ° = 

c 

(c) 

N^hFE"hE^ 

N E 
m 

(d) 

dN  ,,  . . 

Bx  hFE  " hE 

= -f  - (hFE*  F2  layer) 

These  conditions  determine  a cubic  which  has  a minimum  at  the  height 


h 

n 


h E + 
rn 


-h 

m 


E) 


The  complete  profile  is  then: 


(i)  F2  region 


h__  ^hsh  F2 
FE  m 


N (h)  = N F2  J 1 - h F2  -h 
m ( I m 


y F2 
m 


(ii) 


hmE  ‘ h < hFE 


N(h)  = N E + A(h  - h E) 
m m 


(h__  - h E) 
' FE  m 


h - h E 
m 

h_„-h  E 
FE  m 


where  h_„  = h F2  - Jz  y F2  / 1 - (foE  \ 

FE  m V VfbF2 ) 


1/2 


and  A = 


M(*FE>  = iNmF2jl- 
dh 


Jz  y F2 


1/2 


foE 
foF2/  J 


m 


(iii)  E region  h E - >Jz  y E ^ h < h E 


m 


N(h)  = N E 1 + 2 
m 


/ h E - h\ 


m 


\^vE  I 


-3  / h E - h 


m 


Jz  y E 
m 


The  electron  density  in  region  (ii)  at  h is 

n 


N (h  ) 
n 


N E 
m 


77  N FZ[ 
27  m 


h,  -h  E 
FE  m 


y F2 
m 


G 


where 


We  made  the  additional  requirement  that  the  minimum  density,  Nfh^) 

N(h  ) ^ — N E.  If  this  is  not  the  case,  we  set  N(h  ) = — N E and  use 
' n 9 m ' n 9 m 

the  profile  in  this  region  to  redefine  the  semithickness  of  the  F2  region 
obtaining 
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1 


y F2  = - 

m Ji 


N F2  . G . (h  F2  - h E) 
m m m 

3 N E + G(i) 2  N F2 
m m 


The  final  case  is  when  all  three  layers,  F2,  FI  and  E are  to  be 
present.  Here,  we  also  require  that  the  electron  density  within  each 
layer  be  consistent  with  that  predicted  independently  for  the  contiguous 
layer.  The  thickness  of  each  layer  is  defined: 


y = h F2  - h FI 
cm  m 


and 


Y 


D 


h FI 
m 


h E 
m 


A three  layer  profile  is  used  when  y 2:  y F2.  If  this  condition  is  not 

c m 

satisfied,  we  assume  that  the  electron  density  in  the  F2  region  cannot 
fall  fast  enough  to  account  for  the  existence  of  an  FI  layer  at  its  pre- 
dicted position.  When  this  occurs  and  foFl  ^ foE,  we  neglect  the  FI 
layer  and  use  a two  layer  monotonic  profile  as  was  previously  described. 
If  foFl  < foE,  we  still  omit  the  FI  layer  but  use  the  nonmonotonic  profile 
of  an  E and  F2  layer. 

The  three  layer  profile  consists  of  the  following: 


(i)  F2  region  h FI  ^ h s h F2 

mm 


N (h) 


N F2  + (N  F2-N  FI) 
m mm 


(ii)  FI  region 


h ESh<h  FI 
m m 


N(h) 


N FI  + (N  FI  - N E) 
m mm 


2 h FI  - h 
m 

yD 


\3 


h FI  - h 
m 

l y D 


(iii)  E region  h E - Jz  y Es  h <h  E 

m m m 


N(h)  = N E 
m 


+ 2 


/ h E -h] 
m 

y E 
' m 


3 


' h E - h\ 
m 

Jz  y E , 

' m ’ 


Thesa  are  three  cubic  profiles  with  = 0 at  each  layer  height  of 
maximum  electron  density.  No  other  conditions  are  placed  upon  the 
layer  parameters  other  than  those  already  stated  so  that  the  profile 
need  not  be  monotonic.  It  is  possible  that  the  maximum  electron 
density  occur  at  either  the  FI  or  E layer  maximum  height.  A sample 
profile  is  shown  below. 
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CAL  ELECTRON  DENSITY  PROFILES 
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Figure  III- 1 


Figure  III- 1 illustrates  two  vertical  electron  density  distributions  for 
the  case  in  which  foFl^O.  The  parameters  used  in  the  calculation  of  the 
electron  density  at  10  km  intervals  are  the  following: 


Parameter 

Solid  Curve 

Dashed  Curve 

Eocal  Time 

16.  7 

13.  3 

Latitude 

26.  9° 

31.3° 

Longitude 

341. 2° 

289.  1° 

Sec  y 

-5.  87 

6. 42 

fo  E 

3.  0 MHz 

4.  65  MHz 

foF2 

5.  2 

8.  7 

y F2 
m 

47 . 6 km 

85.4  km 

h E 
m 

113.6 

117.  9 

h F2 
m 

312.  7 

289.  0 

h 

min 

194.  0 

163.  9 

The  value  N(h  ) for  the  solid  curve  is  constrained  to  be  1/9  N E. 

All  XI  1 

In  many  applications,  the  vertical  electron  density  model  is  used  to 
generate  a three-dimensional  electron  density  distribution  in  space  and 
time.  This  requires  that  the  longitudinal  gradients  generated  by  using 
the  model  at  different  points  in  space  be  consistent  with  those  that  are 
observed.  Since  we  may  use  an  entirely  different  profile  depending 
upon  the  values  of  ionospheric  parameters,  care  must  be  used  in 
specifying  the  closeness  of  the  calculated  profiles.  Examples  of  two 
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dimensional  ionospheres  is  shown  in  Figures  HI-2  and  III- 3 where  contours 
of  constant  electron  density  are  shown  as  a function  of  height  above  the 
earth  and  angular  range  for  a great  circle  path  of  50  . In  both  figures,  the 
electron  density  is  calculated  at  height  intervals  of  10  km  using  our  vertical 
electron  density  distributions.  In  Figure  ILL- 2,  the  profiles  are  calculated 
at  angular  intervals  of  5°  whereas  in  Figure  III- 3 the  interval  is  1 . The 
isoionic  contours  are  then  drawn  using  a quadratic  smoothing  technique 
using  the  respective  grids  of  (11.  33)  and  (51,  33)  data  points. 

Although  most  of  the  features  contained  in  each  of  the  figures  are 
similar,  there  are  some  regions  in  which  large  differences  exist.  The  mos. 
notable  of  these  are: 

a)  Lower  E region,  9 = 4,  height  = 150  km. 

b)  Lower  E region,  9 = 18  , height  = 125  km. 

o 

c)  Lower  F region,  9 = 38  , height  = 200  km. 

The  gradients  in  each  of  these  regions  for  the  ionosphere  of  Figure  III- 3 
is  a factor  five  times  that  of  Figure  III- 2 . This  is  due  to  the  interpolation 
of  contours  of  constant  electron  density  from  vertical  profiles  which  are 
specified  at  intervals  that  are  five  times  more  closely  spaced.  If  the  interval 
were  more  closely  spaced,  these  gradients  would  become  ever  larger.  This 
deficiency  in  the  model  can  be  surmounted  by  matching  the  fintnes'-  of  the 
choice  of  the  grid  spacing  with  the  magnitude  of  the  gradients  of  electron 
density  occurring  in  the  ionosphere. 
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CONTOURS  OF  CONSTANT  ELECTRON  DENSITY 
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IV.  COORDINATE  CONVERSION  TECHNIQUE  FOR  OTH  BACKSCATTER  RADAR 

Detailed  documentation  of  the  program  WIMP  and  its  ionospheric  and 
mathematical  techniques  has  been  prepared  for  the  Electromagnetic 
Sciences  Division  of  the  Deputy  for  Electronic  Technology.  This  writeup 
along  with  an  annotated  program  listing  is  contained  in  a separate  report. 


The  WIMP  Program  is  a complex  algorithm  which  performs  an  accurate 
conversion  from  radar  target  range  coordinates  to  geographic  target  co- 
ordinates in  OTH  Backscatter  HF  radar  applications.  The  complexity  is 

I 

due  to  the  use  of  a iterative  algorithm  that  involves  interaction  between 
the  program  and  the  user.  The  program  consists  of  three  main  component  . 


a.  A set  of  subprograms  which  construct  a three-dimensional 
three-layer  ionospheric  model. 

b.  A set  of  subprograms  which  simulate  HF  radio  wave  propagation 
in  the  ionosphere  constructed  in  (a). 


c.  One  of  three  driving  programs  controlling  the  calling  sequence 
to  the  radio  wave  propagation  simulator  in  (b).  These  driving 
programs  perform  the  following  tasks:  (1)  From  a given  iono- 


sphere, generate  the  leading  edge  of  a simulated  oblique  ionogram 
including  scale  factors  to  the  F^  layer  parameters  to  lorce  agree- 
ment with  a given  vertical  ionogram;  (2)  from  a comparison  of  the 
simulation  in  (1)  to  a given  oblique  ionogram  generate  range 
gradient  factors  to  apply  to  f F and  M(3000)F  to  force  agree- 

O b 6 

ment;  (3)  from  the  final  ionosphere  generated  in  (1)  and  (2), 


simulate  predicted  radio  frequency  propagation  paths,  from 
which  ranee  versus  eroup  path  functions  mav  be  tabulated. 


The  writeup  contains  a description  of  the  WIMP  three-dimensional 
three-layer  ionospheric  model,  followed  by  a general  description 
of  the  radio  wave  propagation  simulation  subprogram  (ray-tracing 
algorithm).  A functional  description  is  given  of  the  three  driving 
programs  and  the  operating  instructions  for  their  use  are  included. 

Program  block  diagrams,  flow  charts,  and  program  listings  are 
included  in  appendixes. 

Comparisons  have  been  made  of  the  results  of  the  WIMP  ray 
tracings  to  techniques  developed  at  RADC.  The  description  and 
results  of  these  studies  are  reported  in  a separate  annex  to  this 
report. 
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RADC  plans  and  conducts  research , exploratory  and  advanced 
development  programs  in  command,  control,  and  coemunlcatlons 
<C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
ara  communications , electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microeeve 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 
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